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We theoretically analyze the synchronized and desynchronized phases of coupled non-equilibrium 
polariton condensates within mean field theory. An analytical condition for the existence of a 
synchronized phase is derived for two coupled wells. The case of many wells in a 2D disordered 
geometry is studied numerically. The relation to recent experiments on polariton condensation in 
CdTe microcavities is discussed. 
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Polariton condensates in semiconductor microcavities 
provide us with novel macroscopic quantum objects, 
whose precise experimental [l[ and theoretical 0, 0, 
[B|, [|| characterization is presently an active field of re- 
search. A realistic description of the actually realized 
polariton condensates requires to take into account the 
interactions between polaritons, pumping, losses and the 
disordered potential landscape. This potential is due to 
fluctuations in the growth process of the semiconductor 
heterostructure and contains in CdTe microcavities typ- 
ically a few deep minima within the condensate area. 
In the presence of such a potential, the question natu- 
rally arises whether a single condensate is formed that is 
modulated by the disorder potential yet phase-coherent 
over its size, or rather several independent condensates 
with different frequencies are formed. In recent experi- 
ments [2j, it was found that depending on the configura- 
tion of the external potential (position on the sample), 
both the case of a fully coherent condensate with a sin- 
gle frequency (synchronized) and the case of incoherent 
condensates with different frequencies (desynchronized) 
can be realized. The experimental data suggest that the 
underlying physics is analogous to the technologically im- 
jortant phenomenon of mode locking in laser gyroscopes 
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In the present paper, we will present a theoretical pic- 
ture of the mode synchronization in polariton conden- 
sates. Because the average frequency of the polariton 
condensate is related to the macroscopically occupied 
state, a mean field description can be used in a first ap- 
proximation. A non-equilibrium mean field model for po- 
lariton condensates was proposed in Ref. [2| and used to 
describe their peculiar spatial and spectral shape in Ref. 

A similar model was introduced in Ref. (if. Within 
the latter framework, mode locking effects between dis- 
crete energy levels polariton condensates are investigated 
by P. Eastham 0. 

The simplest system that allows to understand the syn- 
chronization physics of a polariton condensate in a disor- 
dered microcavity consists of two coupled wells with a po- 
tential difference. Our analysis will show that a Joseph- 



son flow is responsible for the phase locking between the 
condensates. If the potential difference between the wells 
exceeds a certain critical value, the Josephson flow cannot 
reach a steady state and the condensates no longer share 
the same frequency. Two condensates with different fre- 
quencies are formed in each well and density oscillations 
reminiscent of the AC Josephson effect appear. 

The dynamics of the condensate macroscopic wave- 
function tp{ r ) Wln be described by a generalized Gross- 
Pitaevskii equation [2| ■ Following the work on Josephson 
physics with spatially se par ated ultracold atomic Bosc- 
Einstein condensates [ll|, 13 1 , the condensate wave func- 
tion is projected on the wave functions 4>i^ of each well 
(normalized as usual as Jdr \<f>j | 2 = 1). In terms of the 
amplitudes ipi t 2 in the two wells, the total polariton wave- 
function reads ^r) = tpi (j)i(r)+t/j2 <^2(r) and the dynam- 
ics is given by [2J 



dt 



.[i?(n,)- 7 ]V. 
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(1) 



where 7 is the polariton line width (for simplicity taken 
to be equal in both wells), J is the hopping energy and tj 
is the ground state energy of each well in the absence of 
coupling. The term R(rij) describes the gain of the con- 
densate due to the stimulated scattering from the exci- 
tonic reservoir into the lower polariton states. Polariton- 
polariton interactions are described by the charging en- 
ergy Uj\ipj\ 2 and the interactions between the reservoir 
and condensate polaritons are taken into account by the 
term U^rij. 

The reservoir occupation ni,2 results from the equilib- 
rium between the pumping at a rate Pi 2 and its decay: 



drij 
~df 



Pj -jRn 3 - R(nj)\ipj\' 



(2) 



where the term JrUj represents losses through channels 
different from the condensate. The stationary state and 
excitation spectrum of the symmetric system (ei = £2, 
XJ\ = TJ% and Uf 1 — U^) are discussed in Ref. [2], where 
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it was found that the Josephson plasma oscillations are 
damped due to the non-equilibrium nature of the polari- 
ton condensate. 

Let us start to analyze under which conditions a syn- 
chronized steady state exists if the two potential wells 
arc detuned. The wave functions in the two wells then 



take the form ^1,2 (*) 



-wt±4 



' and rijlt) = rij 



Substituting this state in the equations |T]) and ([2]) gives 

^VpJ = ( e i + u iPj + u P n i)VPj 

-exp[-i(-l 
Pj = iR^j + R(nj)Pj 



J^plZ]exp[-i(-iye} + -[R( nj ) - j}^p] (3) 

(4) 



The imaginary part of Eq. ([3]), expressing the conserva- 
tion of polariton density 



[7 - i?(ni i2 )]pi,2 = ±2J^fp 1 p 2 sin(0), 



(5) 



shows that a Josephson current Ij = IJjpxpi sin(0) 
flows between the two condensates. With the help of 
Eqns. ((3]) and ([S]), the recombination rate R(ni y2 ) can 
be written as 



R{ni. 2 ) = 7 - [Aw tan 9 T y (Aw tan 9) 2 + 4 J 2 sin 2 9 ] , 

(6) 

where the effective detuning Aw is shifted from the bare 
detuning Ae = e\ — e 2 by the mean field shifts: 

Aw = Ae + Utp 1 + U?m -U 2 p2-U?m. (7) 

Eq. © together with Eq. (0| finally allows to find the 
phase difference 9 as a function of the effective condensate 
detuning Aw. 

A simple analytical expression for Aw as a function of 
the phase difference 9 exists for P\ = P 2 = P and when 
the decay of reservoir polaritons is most efficient through 
the condensate: R(rij)pj 3> TrTIr, a condition that is 
expected to hold well above the threshold. The solution 
to Eqns. ^ and © then reads 

2J 2 

Aw = sin(20), if R{nj)pj > -y R n R . (8) 

7 

As expected, the maximal value of the detuning Aw c = 
2 J 2 /7 increases as a function of the coupling between the 
wells J and is inversely proportional to the polariton line 
width 7. 

Under the condition that the nonlinearity Uj in both 
wells is equal, also a simple analytical expression for the 
bare detuning Ae [see Eq.©] as a function of the phase 
difference can be derived. Some algebra with Eqns. © 
to © yields 



2J 2 . (om UU lP ° 
sin(20) 

7 7 
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1 + (2 J/7) 2 sin 2 9 
U^ ni (9) + U^n 2 (9), (9) 



where p° is the condensate density in both wells in the 
absence of coupling (J = 0) and ni :2 (9) is defined by 
Eqs. © and ©. Note that for vanishing tunneling rate 
J, a synchronized solution only exists for Ae = 0: local 
interactions alone cannot lock the two spatially separated 
condensates to a single frequency. 

A plot of the detuning as a function of the phase dif- 
ference between the two wells according to Eq. © is 
shown in Fig. Q] (a) for several values of U\p° and zero 
interaction strength between the reservoir and conden- 
sate U^ 2 = 0. Up to a critical value of the detuning, a 
stationary state with a single frequency for the two con- 
densates exists. The detuning reaches a maximum value 
Ae max for a phase difference 9 w 7r/4, at which the flow 
I j = 2 J^/ P ip 2 sin(#)is maximal. For a detuning larger 
than Ae max no stationary synchronized solution exists. 
It is important to note that the stationary synchronized 
state does not coincide with a linear eigenstate of the 
two- well system, for which 9 S {0, 7r}, but that the con- 
densation occurs in a new state that is formed above the 
threshold for condensation. 

The existence of a synchronized state up to a critical 
detuning is very similar to the phenomenon of mode lock- 
ing in lasers that is most simply described by the Adlcr 
equation for the phase difference 9 between two modes 
[8J. The most important difference between polariton 
condensates and ordinary lasers is the large value of the 
polariton nonlinearity as compared to the very small pho- 
ton nonlinearity in ordinary lasers. Fig. [T](a) shows that 
the polariton interactions help the synchronization of the 
polariton condensate. The physical reason is that the 
flow depletes the high energy well. As a consequence, 
the blue shift of the higher well decreases and the energy 
levels are pulled toward each other. The phase diagram 
as a function of the detuning Ae and pump power (ex- 
pressed in terms of p°) is shown in the inset of Fig. Q] 
(a). 

The interactions between the reservoir and condensate 
polaritons can enhance the synchronization as well [see 
Fig. [U (b)]. Where the condensate-reservoir interac- 
tions counteract synchronization at small phase differ- 
ence, a strong enhancement is obtained for 9 close to 
it/2. Within the approximation that — 0, the contri- 
bution to the detuning from condensate-reservoir interac- 
tions does not depend on the pump intensity. The exper- 
imental study of the synchronization transition as a func- 
tion of the pump power could therefore give an indication 
whether rather the condensate-condensate or condensate- 
reservoir interactions are the dominant mechanism for 
the synchronization. Unfortunately, the independence of 
the pump intensity is spoiled when the finite line width 
of the reservoir excitons jr is taken into account. An 
alternative physical quantity that could give an indica- 
tion about the dominant mechanism for the frequency 
synchronization of the condensates is the relative phase 
between the two condensates. A numerical stability anal- 
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FIG. 1: (Color Online) (a) The effective detuning Ae (in units of the tunneling rate J) as a function of the phase difference 
9. The damping rates are taken 7/ J = 7_r/J = 1 and equal pump rates Pi = P2. The blue shifts in the absence of the 
coupling Up" are shown in the legend and [7^ 2 = 0. The inset shows the maximum detuning as a function of blue shift in the 
absence of coupling, (b) The same for zero condensate-condensate interaction f/1,2 = and non- vanishing reservoir-condensate 
interactions U R = U R = Ur (see legend). The inset shows the maximum detuning as a function of the blue shift in the absence 
of coupling UrUr. 



ysis of the synchronized state and a full integration of 
the motion equations and |2]) have shown that the 
expected window for the relative phase when condensate- 
condensate interactions prevail is < \9\ < 7r/4, where 
tt/2 < \8\ < 7T when condensate-reservoir interactions 
dominate. 



(a) 

1 \A» — 



4 

CD 2 

« 
=T -2 



5 10 
tJ 



(b) 





w/J 



°o. 0.6 
t£T 0.4 
0.2 




(c) 




V \> 


\ 

\ 


) 1 2 
t J 


3 <■ 


(d) J 








> 2 4 


6 S 



FIG. 2: (Color Online) Time evolution of the density (full (pi) 
and dashed (p2) lines) and phase difference (dash-dotted line) 
in a two-well geometry with frequency detuning Ae/J = 2 (a) 
and Ae/J — 3.8 (c). (b) and (d) show the temporal Fourier 
transforms. Other parameters: 7/ J = 1, jr/J = 5 and 
UP/ J 2 = I. 

In agreement with the analytical prediction, the nu- 
merical integration of Eqns. (pj and ((2|) with the param- 
eters for Fig. [5^,,b shows a synchronized solution. The 
densities evolve to a constant value that is symmetric 



around the density in the absence of coupling p = p° . The 
condensate at higher energy (full line) is less occupied 
due to the flow into the lower energy one (dashed line). 
The temporal Fourier transform (Fig. ^Bp) is peaked at a 
single frequency for both condensates. 

In Fig. Ht,d, the ratio of the frequency detuning with 
respect to the coupling Ae/J = 3.8 is too large to allow 
for the locking of the phases in the two wells. The time 
dependent phase difference induces an oscillating Joseph- 
son current that bears a striking analogy with the AC 
Josephson effect, where the application of a constant volt- 
age (chemical potential difference) also leads to an alter- 
nating current. A crucial difference with the AC Joseph- 
son effect in atomic Bose-Einstein condensates concerns 
the relaxation: where the AC Josephson oscillations of 
an atomic condensate are damped at any finite tempera- 
ture and its steady state is at thermodynamical equilib- 
rium with a single chemical potential, no relaxation to 
a single frequency state is possible for a desynchronized 
non-equilibrium polariton condensate. 

In the frequency domain, the alternating Josephson 
currents cause each condensate to have a small contribu- 
tion at the frequency of the other condensate (see Fig. 
[2}i). These overlapping frequency components give rise 
to residual coherence and thus interference fringes. So 
far, this residual coherence was not observed experimen- 
tally Q in the desynchronized regime. To verify our pre- 
diction of residual coherence and density oscillations in 
the desynchronized regime, the experimental investiga- 
tions could be performed in a more controlled way by 
using polariton condensates in mesa microcavities [14j . 

It should also be mentioned that for values of the de- 
tuning, close to the maximal value Ae < Ae c , numerical 
integration of the Josephson equations shows that both 
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the solutions with stationary and oscillating densities can 
be reached, depending on the initial conditions. 
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FIG. 3: (Color Online) The probability to have a synchro- 
nized condensate in a 2D disorder potential with uniformly 
distributed energy levels in the interval [0, Ae] is shown in 
gray scale. Only interactions within the condensate are con- 
sidered {U\2 = 0). The full line shows the boundary of the 
locked solution for two wells that are detuned by an energy 
Ae. The dashed lines show the border of synchronization 
when only the central well in the 2D geometry is detuned by 
an energy +Ae (lower dashed line) or — Ae (upper dashed 
line). From a a 5x5 square of coupled wells, the central 3x3 
part is uniformly excited by the pump laser. 

The generalization of the model equations fl} and © 
to the case of multiple wells with randomly distributed 
energy levels, closer to the disordered reality of CdTe 
microcavities, is straightforward. Fig. [3] shows the prob- 
ability to reach a synchronized state as a function of the 
disorder strength and interactions. The averaging of the 
synchronization over many realizations of the disorder 
potential gives rise to a transition region in the phase 
diagram where the synchronization depends on the ac- 
tual realization of the disordered potential. As a general 
trend, the interactions increase the probability to reach 
a synchronized state as well as the width of the transi- 
tion region. For comparison, the analytically determined 
phase boundary of the two-level system (JSJ) is shown by 
the full line. The dashed lines show the phase boundary 
for a regular 2D array of energy levels where the central 
level is detuned by an energy +Ae (lower dashed line) and 
by an energy — Ae (upper dashed line). These lines show 
that a higher dimensionality favors the synchronization 
for the non-interacting polariton condensate. A physi- 
cal explanation for this fact is that the total Josephson 
flow out of (or into) the central well can reach a higher 
value because it is distributed over more links. Another 
clear feature is that the interactions help the synchroniza- 
tion much more if the central well is negatively detuned. 



Note finally that the presently considered inherently non- 
equilibrium desynchronized phase is very different from 
the thermodynamic equilibrium insulating phase studied 
for microcavity polaritons in Ref. [f| . 

In summary, we have analyzed the recently observed 
synchronization-desynchronization transition of polari- 
ton condensates within mean field theory. An analytical 
condition for synchronization is derived for the case of a 
two-well system. The many well configuration was ana- 
lyzed numerically. The same phenomenology was found 
in both cases: for small detuning between the different 
wells, the Josephson currents and densities reach a steady 
state. For too large detuning on the other hand, a steady 
state no longer exists and the Josephson currents cause 
density oscillations in the different condensates. Both 
interactions of the condensate with itself and with the 
reservoir are shown to enhance the synchronization. The 
similarities and differences with the Josephson effect of 
superfluids and the locking-delocking transition in ordi- 
nary lasers was clarified. 
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